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ABSTRACT 

We study the effects of optically thin radiative cooling on the structure of radia- 
tively inefficient accretion flows (RIAFs). The flow structure is geometrically thick, and 
independent of the gas density and cooling, if the cooling time is longer than the viscous 
timescale (i.e., i coo i > i V isc)- For higher densities, the gas can cool before it can accrete 
and forms the standard geometrically thin, optically thick Shakura-Sunyaev disk. For 
usual cooling processes (such as bremsstrahlung) , we expect an inner hot flow and an 
outer thin disk. For a short cooling time the accretion flow separates into two phases: 
a radiatively inefficient hot coronal phase and a cold thin disk. We argue that there 
is an upper limit on the density of the hot corona corresponding to a critical value of 
tcooiAff(~ 10 — 100), the ratio of the cooling time and the free-fall time. Based on our 
simulations, we have developed a model for observed transients in black hole X-ray bina- 
ries (XRBs) . An XRB in a quiescent hot RIAF state can transition to a cold black-body 
dominated state because of an increase in the mass accretion rate. The transition from 
a thin disk to a RIAF happens because of mass exhaustion due to accretion; the transi- 
tion happens when the cooling time becomes longer than the viscous time at inner radii. 
Since the viscous timescale for a geometrically thin disk is quite long, the high-soft state 
is expected to be long-lived. The different timescales in black hole transients correspond 
to different physical processes such as viscous evolution, cooling, and free-fall. Our model 
captures the overall features of observed state transitions in XRBs. 

Key words: accretion, accretion disks; X-rays: binaries 



X 

S-H 



1 INTRODUCTION 

Black holes (BHs) are the simplest astrophysical objects, com- 
pletely specified by their mass and spin. Black holes are also 
the most powerful sources in the universe because of their 
strong gravity. Roughly ~10% of the rest mass energy of the 
matter falling into BHs can be extracted in the form of ra- 
diation and jets/outflows. While BHs themselves are quite 
simple, accretion flows around them are not. The accretion 
flows are affected by complex, nonlinear processes such as 
turbulent angular momentum transport, cooling, and radia- 
tion (see Balbus fc Hawley |199 8; Nara yan fe Quataert |2005 
for a review). 
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There are two well-established classes of astrophysical 
black holes: stellar- mass black holes (~ 10 M0), the end 
products in the evolution of some of the massive stars (see 
Woosley, Heger, fe Weaver |2002 for a review); and supermas- 
sive black holes (~ 10° — 10 lu Mq) at the centers of galax- 
ies. Both these classes of BHs show two dominant states of 
the accretion flow around them: an optically thin, geometri- 
cally thick, hot radiatively inefficient state (e.g., see Ichimaru 
|1977 Rees et al. |1982 1 ; and an optically thick, geometrically 
thin black-body disk state (Shakura fc Sunyaev ||1973|. The 



quasars observed at high redshifts, with a dominant black- 
body component, are examples of thin disk accretion in su- 
permassive BHs. The quiescent active galactic nuclei (AGNs) 
at low redshifts lack thermal disk emission and are modeled 
as radiatively inefficient accretion flows (RIAFs). Because of 
shorter timescales in X-ray binaries (XRBs) the same sys- 
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tern can be observed to transition from a quiescent state to a 



black-body state and vice- versa (see Remillard & McClintock 
2006 for a review). 



Most accretion flow simulations model the RIAF state 
because one does not need to worry about cooling and ra- 



diation transport (e.g., Stone et al. 



1999 



Igumenshchcv, 



Narayan, & Abramowicz 2003; Bcckwith, Hawlcy, & Kro 
lik ||2008| Penna et al. ||2010 1. Another problem in modeling 
thin disks is the enormous resolution required for a realistic 
H/R ~ 0.001-0.01 (ratio of disk height to radius). There are 
a handful of simulation papers which include some effects of 
cooling on RIAFs (e.g., Li, Ostriker, fe Sunyaev |2012 Barai, 



Proga, fc Nagamine |2012 |. Progress has been made recently 
to include most of the relevant physics (magnetic fields, radi- 
ation, cooling, etc.) in numerical simulations (e.g., |Ohsuga et| 
al. |2009 Ohsuga fc Mineshige |2011 1 but we are still far away 
from modeling all the processes with realistic parameters. 

In this paper we perform detailed hydrodynamic simu- 
lations of accretion with varying viscosity prescriptions. The 
main aim is to study the influence of optically thin radiative 
cooling on the global structure of the accretion flow. While it 
is widely believed that magnetohydrodynamic (MHD) turbu- 
lence associated with the magnetorotational instability is re- 
sponsible for angular momentum transport and dissipation in 
accretion flows ( Balbus fc Hawley |1991 1 , a phenomenological 
treatment of viscosity allows us to gloss over the complexities 
(both numerical and physical) of MHD accretion and to focus 
on important physical processes such as the viscous accretion, 
cooling, and response to gravity. 

We confirm the theoretical prediction that a thin disk 
forms at radii where the cooling time is shorter than the vis- 
cous accretion time ( Rees et al. |1982 l. We observe an upper 
limit on the density of the hot corona corresponding to a lower 
limit on the ratio of the cooling time and the free-fall time, 
icooi/iff ~ 10 — 100, irrespective of the form and magnitude 
of viscosity. This means that there is an upper limit on the 
amount of mass that can go into the hot corona; the rest of 
the mass has to cool onto the thin disk. 

Motivated by our idealized simulations, we identify dif- 
ferent physical processes responsible for the observed char- 
acteristics of black-hole binary transients. We argue that the 
observed transition of the quiescent (very low-hard) state to 
a black-body dominated (high-soft) state via a short-lived 
high/hard state occurs because of a large increase in the mass 
accretion rate. This transition from hard to soft state can oc- 
cur at a fast timescale, comparable to the viscous timescale 
at the inner circularization radius. The transition from a thin 
disk back to the low-hard state happens only after sufficient 
mass has been accreted, so that the cooling time is longer 
than the viscous time everywhere. This occurs at the longest 
viscous time in the thin disk, which can be years. 

The paper is organized as follows. In section 2 we present 
the hydrodynamic model and details of numerical implemen- 
tation. In section 3 we discuss our results of runs with and 
without cooling. In section 4 we present a model for global 
properties of BH XRB transients based on our simulations. 
In section 5 we conclude with astrophysical implications. 



2 METHOD 

2.1 Basic Hydrodynamic Equations 

In our numerical simulations we solve the basic hydrodynamic 
equations in spherical polar coordinates (r, 6, </>) using the 
ZEUS-MP code ( |Hayes et al. |2006[ ). These are, the equation 
of continuity, 

f t+p V-v = 0, (1) 



the momentum balance equation 



p— = -VP-pV(f> + V -a, 



and the energy equation, 



p dy£ = _ PV . v + a ? /li . 



n e n;A(T). 



(2) 



(3) 



Here, p is the mass density, P the pressure, v = (v r ,ve, v^) the 
velocity, e the internal energy density (= 3P/2 corresponding 
to an adiabatic index 7 = 5/3; protons are non-relativistic till 
close to the BH), p = vp the coefficient of viscosity (y is the 
kinematic viscosity), a the viscous stress tensor, n e and m 
the electron and ion number densities, and A(T) the cooling 
function with temperature T. The Lagrangian time derivative 
is denoted by d/dt = d/dt + v ■ V. 

We use the pseudo-Newtonian potential proposed by 
Paczyhski and Wiita (1980), which describes the gravitational 
field of the central black hole, 



GM 
~-R a 



(4) 



where M is the the black hole mass, c the speed of light, G the 
gravitational constant, and R g (— 2GM/c 2 ) the Schwarzschild 
radius. 

We include a bremsstrahlung cooling term in the energy 
equation (last term on the right hand side of Eq. p}| . We use 
the cooling function A(T") given by equation ( 12) of|Sharma 
2010 1, which has been adopted from ISutherland & 



et al. 



Dopita I ( 1993 1 for mean molecular weights corresponding to 
p = 0.62 and p e = 1.18. We note that in the present work we 
ignore cooling due to synchrotron and inverse Compton pro- 
cesses, which although important in certain accretion regimes 
(e.g., Rajesh fc Mukhopadhyay |2010 |, are difficult to model. 
Both synchrotron and Compton cooling terms are expected 
to be proportional to n 2 (square of number density) because 
both the magnetic and radiation energy densities are expected 
to scale with the mass accretion rate. Thus, including these 
should give results qualitatively similar to bremsstrahlung 
cooling. 

The second term on the right hand side of Eq. H 
represents heating due to viscous dissipation. MHD simula- 
tions (e.g., Table 2 in Hawley et al. |1995 Brandenburg et aL] 
||1995 1 show that the r<f) component of the magnetic stress is 



the most dominant component of the stress tensor. Although 
we do not include magnetic fields in this work, in order to 
mimic its effects, we assume that only the r<f> component of 
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the viscous stress tensor a is non-zero, such that: 



onto the black hole, is given by 



cr r $ = a$ r = M r 



d (v^ 
Or 



(5) 



This anisotropic stress behaves very differently from an 
isotropic viscous stress, and only damps radial variations in 
the angular velocity. 

We choose the fiducial form of the kinematic viscosity 
in our runs to be v oc r 1 ' 2 , guided by the existing theoret- 
ical predictions. In the standard 'a— disk' model of a thin 
disk, v — ac s H, where a is the Shakura-Sunyaev viscosity 
parameter, c 3 the sound speed and H the disk scale height 
( Shakura fe Sunyaev 1 1973 1. Now, H ~ c s /Q, where fi is the 
angular velocity, and we have, v w aH 2 Q fa a(H/r) 2 Q.r 2 . 
For a Keplerian flow, Q oc r~ 3 ' 2 and hence, v scales with r 



as v oc r 



1/2 



provided H/r is constant, which is indeed the 
case for a Shakura-Sunyaev disk. For comparison we also try 
v — constant and v oc p. The kinematic viscosity for the 
fiducial case is given by 



i/ = 4.5x 10" 4 a .oin i?o 



Ro 



1/2 



(6) 



where ao.oi = a/0.01, Ro is the center of the torus where the 
initial density maximum is located (see section 2.3) and ilo is 
the angular velocity at Ro; most runs use a = 0.01 but some 
runs which are evolved for long times use a = 0.1. Therefore, 
the viscous time is sa 16000^^ times the free-fall time (see 
Eqs. U\ & [81. The magnitude of the constant viscosity (for 
v — constant case) is obtained by evaluating v in Eq. 6 at 
Ro- Similarly, the normalization for v oc p case is chosen such 
that the viscous timescale at Ro is the same as in v oc r 



1/2 



Note that, in absence of cooling, Eqs. 1-3 are independent 
of the density normalization (since P oc p and a oc /i oc p, 
provided the kinematic viscosity v is not a function of den- 
sity). Thus, the accretion flow structure in absence of cooling 
is identical if we normalize the density by its initial maximum 
value (this can be seen from Figures 1(e) and 3(a) which show 
density at late times for RIAFs with different initial densi- 
ties). Furthermore, the only spatial scale in Eqs. 1-4 is the 
Schwarzschild radius R g , which is proportional to the black 



hole mass (e.g., Rajesh fc Mukhopadhyay |2010 l. Thus, if we 
scale all the spatial scales with R g , velocities with c, time with 
R g /c, density with the initial maximum value po, P and e by 
poc 2 -, v by cR g , then in absence of cooling (i.e., for RIAFs) 
there is no explicit dependence on the black hole mass in Eqs. 
1-4. Therefore, our results can be scaled to a wide range of sys- 
tems — from accretion flows around stellar mass black holes, 
namely X-ray binaries (XRBs), to those around supermassive 
black holes, namely active galactic nuclei (AGNs). 



2.2 Key Timescales and RIAFs 

Here we list some of the important timescales involved in 
accretion flows with cooling. The viscous timescale of the ac- 
cretion flow, which is the time taken by the matter to accrete 



r 2 

v 



We define the free fall timescale of the flow to be 



tg 



2r 3 
GM' 



The cooling timescale of the flow is given by 



cool 



(-) 



(8) 



(9) 



where e = n e n,iA(T) is the cooling rate. If ks is Boltz- 
mann's constant and n denotes the number density then 
e = 3nfcsT/2 and e = n 2 A(T), and hence £ coo i oc n" 1 . Thus 
the cooling time is shorter for a denser optically thin accretion 
flow. 

If the cooling time is shorter than the viscous time (i.e., 
icooi < ivisc) in RIAFs the accretion flow separates into two 
phases, a coronal phase with long cooling times and a dense 
optically-thick, geometrically thin disk. Both these phases are 
in thermal equilibrium: viscous heating balances black-body 
radiation in the thin disk; viscous heating balances optically 
thin cooling and adiabatic expansion in the coronal phase. 
A dense accretion flow develops these two phases because an 
optically thin flow which cools efficiently is globally unsta- 
ble (e.g., Piran |1978 ). The two-phase corona-disk is globally 
stable. If cooling in the coronal phase dominates heating, ex- 
cess mass can condense onto the thin disk and equilibrium 
is achieved. Similarly if heating in the corona dominates, ex- 
cess energy goes into evaporating matter from the thin disk 
and in launching outflows. Although the corona is in global 
thermal equilibrium, it can become locally thermally unsta- 
ble for short cooling times, and excess mass from the corona 
can condense into cold blobs which fall on the thin disk (see 
3~4land|Sharmaet al. |20J2 |. 



section 



In our simulations we do not include radiative transport 
and hence are unable to simulate a realistic standard thin 
disk. Note that, as the RIAF cools, it becomes dense and ra- 
diation does not simply escape, but diffuses with a short mean 
free path. The gas becomes optically thick before cooling to 
10 4 K (which is the temperature of the stable phase for our 
cooling function; cold thin disk exists at this temperature in 
our simulations). Therefore, in reality, the thin disks should 
be hotter and of lower density compared to those in our sim- 
ulations. We note that our thin disks are severely unresolved 
because photon transport (e.g., Ohsuga fc Mineshige ||201 1 ) 



or magnetic fields (e.g., Machida, Nakamura, & Matsumoto 



|2006[ ) which can provide vertical support for the cooling gas, 
are absent. The present paper focuses on the structure of RI- 
AFs and the conditions for the formation of a thin disk, not 
on its detailed structure. Realistic thin disk simulations with 
radiation transport will be carried out in future. 



2.3 Initial Conditions 

Our initial condition consists of a dense torus having a con- 
stant specific angular momentum in dynamical equilibrium. 



© 2013 RAS, MNRAS 000,[T]fl5 



4 U. Das, P. Sharma 



t=o 



1 orb 



24 orbs 




iog 10 P 
-18 




100 200 300 400 
R/R 



200 
R/R 



400 




400 



P 






-18.5 



-19 



-19.5 



-20 



-20.5 



-21 



21.5 



22 



Figure 1. (a) Initial density contour plot for the run with po = 10 6 m p g cm 3 . Panels (b,c) show evolution without viscosity and panels 
(d,e) with viscosity v <x r 1 ' 2 and a = 0.01 (see Eq. Ifjk. The grid resolution is 512 2 . 



The torus is embedded in a medium having very low density 
and is non-rotating. The pressure and density in the torus 
satisfy the polytropic equation of state P — Kp 1 , K being 
a constant. Then following [Papaloizou & Pringl e~||1984| and 
IStone et al. |[T999l we determine the the initial structure of 
the torus, for our potential <f), as: 



P 



GM 



(n+l)R 



Ro 



R„ 



1 

2p^ 



Ro 



R g )R] 



1 
2rf 



(10) 



where, R = rsin#, is the cylindrical radial coordinate; Ro is 
the equatorial distance of the center of the torus (where it 
has maximum density) from the black hole; d (> 1) is the 
distortion parameter which determines the shape and size of 
the torus — if d « 1, then the torus has a small, nearly circular 
cross-section and as d becomes much greater than 1, the cross- 
section deviates from a circular one and the torus becomes 
fatter; n — 1/(7 — 1) is the polytropic index. For most of our 
runs the distortion parameter is taken to be d — 1.125. The 
density of the torus is maximum (p = po) at r = Ro and 
9 = 90°, using which one can write down the expression for 
the constant K: 



K = 



GM 



(n + 1)-RoaT 



Ro 



1 



R{) 



Ro — Rg 2 (Ro 



1 

2d. 



Thus, knowing po and d one can easily determine K. 

The initial density of the embedding ambient medium is 
chosen to be p am b = 10~ 4 po, which is too small to affect our 
results. The mass of the bl ack hole is taken t o be that of Sgr 
A*, i.e., M = 4 x 10 6 M Q (JGhez et al. |2008|). We scale most 



(If) 



of our results for both XRBs and AGNs since the governing 
equations can be scaled for different BH masses. 



2.4 Numerical Setup 

We carry out two-dimensional, axisymmetric, hydrodynamic 
simulations in spherical (r, 8, <fi) geometry, with the computa- 
tional grid extending from an inner boundary i? m m = 1.75i? 9 
to an outer boundary _R m ax = 400R g . A logarithmic grid is 
used along the radial coordinate, while a uniform grid is used 
along the 6 direction, where ^ 9 ^ n. Although there is no 
variation along the <f> direction, it ranges from Sj 4> Sj 2-7T. 

As for the boundary conditions, we use inflow condition 
at the inner radial boundary, such that mass is allowed to 
flow into the black hole but not out of it. We apply outflow 
condition at the outer radial boundary, such that mass is 
allowed to flow out of the grid but not into it. In addition, 
the stress a rl p is set to zero at the inner radial boundary. 

Our standard resolution is N r — Ng = 512, where N r 
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Table 1. Important runs with cooling 



Runs 


a 


Pol m p (cm 3 ) 


Resolution 


A.01n9 


0.01 




10 9 


512 2 


A.OlnlO 


0.01 




10 10 


512 2 


A.Olnll 


0.01 




10 11 


512 2 


A.ln9 


0.1 




10 9 


256 2 


A.lnlO 


0.1 




10 10 


512 2 


A.lnll 


0.1 




10 11 


512 2 



and Ng denote the number of grid points in the r and 8 di- 
rections, respectively. For some simulations, which are evolved 
for many viscous times at large radii and hence require a rela- 
tively longer computational time, we have used low-resolution 
models with N r = N e = 256 or N r = N e = 128. 

In some of the high resolution runs with cooling the den- 
sity can become unusually low at some isolated grid points. 
In order to avoid this and the associated numerical problems, 
we have used appropriate density and internal energy floors 
in our simulations, which however, do not affect the overall 
accretion flow. 



3 RESULTS 

In addition to runs with cooling which are the focus of this 
paper, we carry out a few runs without cooling to study the 
pure RIAF state and its dependence on various parameters. 
The important runs with cooling are summarized in Table 1. 



3.1 Runs without Cooling 

Here we discuss briefly the pure RIAF runs where cooling is 
switched off (i.e., the last term in Eq. (|3j is dropped). As 
discussed earlier, in this regime the density just scales with 
the initial maximum density, and the flow is self similar, if the 
spatial and temporal scales are normalized by the BH mass. 
For pure RIAF runs the maximum density of the torus 
is chosen to be po = 10 6 m p g cm -3 (m p is the mass of the 
proton) at the torus center Ro = 100i? 9 . The density is loosely 



motivated by Sgr A* observations (Baganoff et al. 2003). The 
torus is evolved for about 24 orbits at r = Ro to ensure that 
a steady state is reached (at least at inner radii). Figure [na) 
shows the density contour plot for the initial condition (t = 
0). Since viscosity is responsible for accretion, in absence of 
viscosity (i.e., for v = 0) the torus should not evolve and 
must retain its initial condition. Figures nTb,c) which show 
the density at late times indicate that this is indeed the case; 
a very low density film forms around the initial torus at late 
times as the equilibrium adjusts slightly. This case shows that 
numerical viscosity is negligible and accretion happens only 
because of explicit viscosity in Eq. [2] 

The run with v oc r ' 2 shows viscous accretion and the 
formation of a geometrically thick RIAF at small radii at 
late times (see Figs. fTfd,e)). The matter tapers into a tail 
like structure in the outer region, which is necessary because 



some matter has to move out to carry the angular momentum 
outwards. The inner bulge is a low density (p ;> 10 -20 g cm~ 3 ) 
corona. We use the term 'corona' quite loosely for the hot 
phase of the accretion flow with a long cooling time. Note 
that this 'corona' is formed as a result of accretion and should 
not be confused with the very low density ambient medium 
(Pamb = 10 _4 po), which is required for numerical reasons. 

In order to investigate how sensitive the accretion flow 
is to different parameters, we have carried out several sim- 
ulations by taking different initial values of the distortion 
parameter d, the location of the torus center Ro, and dif- 
ferent functional forms of v. First, the initial torus was 
placed at different locations from the black hole, namely at 
Ro = 50, 100 and 200i? 9 (note that d = 1.125 and v oc r 1/2 
for all three cases). For all the cases, the density contours in 
the inner region of the accretion flow in steady state, were 
found to be very similar in overall shape. Next, different ini- 
tial cross-sectional areas of the torus were chosen, namely 
d = 1.05, which correspond to a smaller nearly circular torus, 
d = 1.3, which consists of a fatter torus deviating from a 
circular cross section, and d = 1.125, which correspond to 
a torus midway between the two extreme cases (note that 
-Ro = 100-R g and v oc r 1 ' 2 for all three cases). Again all the 
cases were found to exhibit very similar density contours in 
steady state. These results indicate that the final steady state 
is not too sensitive to the initial conditions and we conclude 
that the accretion flow is independent of the initial distribu- 
tion of matter around the black hole, and the way the black 
hole is fed, as long as the circularization radius (Ro) is far 
away from the Schwarzschild radius. 

We have also used different functional forms of the kine- 
matic viscosity, namely, v = constant and v oc p (note that 
d — 1.125 and Ro = 100R g ). In order to assess the depen- 
dence of the accretion flow on v, one can arrive at a rough 
scaling of v and density with radius from the conservation of 
angular momentum (c. f. Eq. 2.8 of Pringle |1981 1. In steady 
state this implies that the viscous and the advective fluxes in 
the equatorial plane satisfy 



vpHR 



dQ, MQ.R 1 



(IR 



2k 



(12) 



where the accretion rate M = constant. Now, on assuming 
that the rotational velocity is nearly Keplerian, the above 
equation leads to vpH = constant. Thus the structure of a 
RIAF depends on the form of the kinematic viscosity. The 
results from these two runs (v =constant and v oc p) look 
very similar, as explained in Stone et al. (19991. They differ 
from the v oc r 1 ' 2 case in that, they show a more pronounced 
corona at inner radii. Our results with varying viscosity (v) 
are consistent with Eq. |12| 



3.2 Runs with Cooling 

In this subsection we discuss in detail the effects of 
bremsstrahlung cooling on the accretion flow. Simulations 
without cooling are invariant to any change in po, but the 
runs with cooling are not. Low density runs are essentially 
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Figure 2. (a) (b) (c) (d) The time evolution of temperature contour plots for the fiducial cooling run with po = 10 10 m p g cm -3 and 
a = 0.01. (e) The variation of the mid-plane aspect ratio (H/R) with radius (R/R g ). The black (dot-dashed), blue (dotted), red (dashed), 
and green (solid) lines correspond to the snapshots (a), (b), (c), and (d) respectively. The grid resolution is 512 2 . 



RIAFs throughout (even with cooling), but an increase in 
density eventually leads to the condensation of the hot gas 
and the formation of a cold thin disk. Ideally the accretion 
rate M (or even better M/MecM, where MEdd is the Edding- 
ton rate) should be chosen as the parameter to study this 
transition, but po is a useful proxy. 



3.2.1 The Fiducial Cooling Run 

Figures 2(a)-(d) show the variation of the two-dimensional 
temperature contour plots with time for the accretion flow 
having an initial maximum density po = 10 10 m p g cm -3 and 
viscosity parameter a = 0.01 (see Eq. p] Run A.OlnlO in Ta- 
ble 1). The accretion flow starts off as a geometrically thick, 
hot flow. It gradually develops a peanut-like shape after about 
7 orbits (roughly the cooling time at Rq) and also shows the 
onset of cold high density gas in the mid-plane. After 9 or- 
bits the low density corona displays a pinch at about 100-R 9 , 
and a more extended cold, thin disk is formed. Once steady 
state is attained, one sees the clear formation of an inner low 
density, hot (T ~ 10 10-11 K), geometrically thick corona and 



10" 



K) disk which extends 



a geometrically thin, cold (T > 
till very close to the black hole. 

Since the cooling time is inversely proportional to the 
density, a higher density (or equivalently a higher mass ac- 



cretion rate) corresponds to a shorter cooling time and an 
increase in the radiative efficiency of the flow. Thus, we con- 
clude that the initial density of the torus in this case is high 
enough to trigger the onset of bremsstrahlung cooling, which 
eventually condenses the hot flow to form a geometrically 
thin cold disk. Some corona models, unlike ours (e.g., Meyer 



fc Meyer-Hofmeister |198T I , require conductive heating of the 
upper layers of the thin disk, which evaporates mass from the 
disk to the corona. Our corona is maintained by viscous heat- 
ing; in reality, magnetic dissipation and evaporation due to 
thermal conduction (magnetized thermal conduction is quite 
non-trivial to model faithfully) should also play a significant 
role. Irrespective of the source of heating, the corona exists 
in a state of rough dynamical and thermal equilibrium, where 
maximum density is governed by the physics of local thermal 
instability in presence of gravity (see also section 3.4). 

The transition to a thin disk is further illustrated in Fig- 
ure 2(e), which shows the variation of the disk height to ra- 
dius ratio (H/R; a measure of geometrical thickness of the 
flow) as a function of radius (R/R g ). The disk scale height 
H fa c s /n is calculated using the mid-plane (9 — it/2) value of 
c s -jpjp and n fi^GM/R 3 . The black (dot-dashed), blue 
(dotted), red (dashed) and green (solid) lines show the H/R 
of the accretion flow at 3.5, 7, 9 and 24 orbits, respectively. 
Initially the flow is geometrically thick throughout, with a 
high H/R ;> 0.1. As accretion continues, H/R starts dropping 
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Figure 3. (a) (b) (c) Top row - shows the density contour plots in the steady state for runs with three different initial densities, (d) (e) 
(f) Middle row - shows the contour plots of the ratio of cooling to viscous time (t coo i/t v isc)- (g) (h) (i) Bottom row - shows the variation 
of H/R with R/R g . The grid resolution is 512 2 . 



near the density maximum (around 1007? 9 ), leading to the 
formation of a geometrically thin disk. We note that the in- 
nermost region is always geometrically thick with H/R J> 0.1, 
however, as time progresses, the geometrically thin disk with 
H/R ~ 5 x 10~ 4 moves closer in to the black hole (up to 
about 20R g ). Note that Figure 2(e) shows spikes in H/R in 
the region of the cold thin disk because there is hot gas at 
some of the grids in the equatorial plane (at 2 = 0). This 
increases the estimate of H/R locally. Also note that the ex- 
treme outer region is not a thin disk because the low-density 
hot gas has a long cooling time. 

The accretion flow is expected to have an inner RIAF and 
an outer cold disk because the cooling time, t coo i oc T 1 ' 2 /n, 
is expected to decrease outward faster than the viscous time, 
ivisc = r 2 /v. As we discuss in the next sections, this ratio 
determines whether a thin disk forms at a given radius. 



3.2.2 Dependence on Initial Density 

The fact that a minimum threshold density is required to 
trigger efficient bremsstrahlung cooling is better illustrated 
in Figure 3, which compares results from simulations hav- 
ing three different initial densities. Figure 3 shows the steady 



state contour plots of density (top row), the ratio i coo i/ivisc 
(middle row) , and the variation of H/R with radius (bottom 
row), for accretion flows with different initial density nor- 
malizations: po = 10 9 m p g cm -3 (left column; Run A.01n9), 
10 10 m p g cm -3 (middle column; Run A.OlnlO) and 



po 
Po 



10 11 m p g cm 



(right column; Run A.Olnll). Note that 
all the three cases have the same viscosity parameter a (= 
0.01; see Eq.[6|, i.e., the same viscous time, but different cool- 
ing times due to the difference in the density of the accreting 
matter. 

We observe that the accretion flow with the lowest initial 
density po = 10 9 m p g cm -3 in Figure 3(a) is throughout geo- 
metrically thick with H/R ;> 0.1 (Figure 3(g)) and resembles 
the flow without cooling in Figure 1(e). The entire flow can 
be classified as a hot, low density corona which is radiatively 
inefficient. However, as po increases to 10 10 m p g cm -3 (Fig- 
ure3(b)) the hot corona is restricted to the innermost radii 
and there is appearance of a cold, geometrically thin disk 
with H/R ~5x 10~ 4 (Figure 3(h)). The time evolution of 
this run is already discussed previously (see Figure [2J) . For 
po = 10 m p g cm -3 the cold disk extends very close to the 
black hole (up to <C 10i? 9 ) (Figure 3(i)) and the inner coronal 
bulge has diminished significantly relative to the disk (Figure 
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Figure 4. Temperature contour plots for different runs illustrating the importance of the ratio t C oolAvisc : ( a ) a = 0.01, po = 10 % 
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3(c)). Thus, po = 10 10 m p g cm -3 seems like a threshold den- 
sity, which marks the transition from a geometrically thick 
hot flow to a geometrically thin cold disk. 

The condition for the formation of a geometrically thin 
disk is governed by the ratio of the cooling time and the vis- 
cous time (icooi/ivisc). Figure 3(d) shows that t C ooi/tvisc is 
everywhere greater than 1. Thus the flow is dominated by 
viscous effects and cooling does not affect the flow because 
the gas accretes before it can cool. As the initial density 
(po) increases the cooling time becomes shorter and cooling 
starts to become important. Figure 3(d) shows that the ratio 
icooi/ivisc will become smaller than unity close to the dense 
torus if the initial density is increased by a factor of ten. In- 



deed, the run with po 



10 1 



g cm shows a thin disk 



with tcooi/ivisc <C 1 beyond 20R g in the steady state and the 
inner RIAF has i CO oi/ivisc > 1 (Figure 3(e)). A similar con- 
clusion is drawn from Figures 3(c) and (f) for the case with 
po = W 11 m p g cm~ 3 , which is a completely cooling domi- 
nated flow, resulting in a thin accretion disk extending till 
within 10R g . The H/R ratio in the thin disk for the two high 
density cases is similar (~ 5 x 10 -4 ). As mentioned earlier, 
this ratio for the cold thin disk is not physically consistent 
because vertical radiative transport is ignored in our simula- 
tions. 



3.2.3 Dependence on a Parameter 

So far we have only considered flows with different initial 
densities but with the same viscosity parameter a (see Eq. 
™. If the criterion for the formation of a cold thin disk is 
icooi/ivisc < 1 then increasing a (i.e., decreasing i v isc) should 
result in a higher initial density for which the transition of 
a RIAF to a thin disk happens. To test this, we study two 
runs with a = 0.1 (ten times the fiducial value), and initial 
densities po = 10 10 m p g cm -3 and po = 10 11 m p g cm -3 . Since 
the viscous time is ten times shorter than the fiducial case, 
icooi/ivisc in steady state for po = 10 10 m p g cm -3 and a = 0.1 
run is expected to be similar to the run with po = 10 9 m p g 
cm" 3 and a — 0.01 (Fig. 3(d)). Therefore, we do not expect 



the po 



10 1 



g cm run with a = 0.1 to show a thin 



disk in the steady state. This is indeed the case as can be 
seen from Figure 4(b). 

Figure 4 shows the two-dimensional temperature contour 
plots of runs with different viscosity parameters (as) and ini- 
tial density normalizations (pos) . The left column corresponds 
to the runs with a = 0.01, while the right column corresponds 
to those with a — 0.1. Interestingly, although Figures 4(b) 
and (c) correspond to the same po, they look absolutely dif- 



ferent. After 9 orbits the run with po = 10 



g cm and 



a = 0.01 (our fiducial cooling run) shows the condensation 
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of a cold geometrically thin disk with T ~ 10 4 K, while the 
run with the same density but with a = 0.1 shows a pure 
RIAF. Figure 4(d) shows the formation of a cold thin disk for 
a ten times higher density run (po = 10 n m p g cm -3 ) with 
a = 0.1. Although Figures 4(c) and (d) are plotted at 9 and 
0.9 orbits respectively, they correspond to the same viscous 
time and hence the thin disks look quite similar in the two 
cases. Similarly, while the RIAF temperature contour plots 
in Figure 4(a) and 4(b) are plotted at the same time, their 
viscous times differ by a factor of 10. Therefore the higher 
viscosity run has a more extensive corona; we expect a simi- 
lar flow structure for the a — 0.01, po — 10 9 m p g cm" 3 run 
at 90 orbits. 



3.3 Cooling Induced State Transitions 

There are many observational evidences of black hole bi- 
nary systems that exhibit more than one X-ray spectral state 



(Remillard & McClintock 2006 Fender, Belloni, & Gallo 



20041. The spectra from accretion disks surrounding these 



BHs switch from a hard X-ray dominated state to a disk- 
dominated soft state, and back to a hard state on a timescale 
of months to years. Supermassive BHs in galactic centers are 
also expected to undergo such transients but on much longer 
timescales because the basic timescales are proportional to 
the BH mass. The black-body dominated soft state is well-fit 
by the standard Shakura-Sunyaev thin disk and the hard state 
can be described by a hot RIAF. The transition from a RIAF 
to a thin disk state can be induced by adding mass to the 
accretion flow, thereby increasing the density and decreasing 
the cooling time of the accretion flow such that the critical 
ratio icooi/ivisc becomes less than unity. The transition from 
a dense thin disk to a RIAF occurs naturally once sufficient 
mass in the dense accretion flow is accreted and i coo i/t v isc 
becomes ;> 1 in the inner regions of the accretion flow. 



decrease with time because of accretion. Therefore, to pro- 
duce a transition from a RIAF to a thin disk state we need 
to add mass to the accretion flow. To mimic the increase in 
the mass accretion rate we include a source function in the 
density evolution equation (Eq. nT) given by 

(rsin6>-i?o) 2 + (rcos6>) 2N 



p(r,0)=9 



PQ 

Cnrb 



exp 



(wR g y 



(13) 



where, po is the initial density normalization, and _Ro is the ra- 
dius where mass addition is centered (chosen to be the center 
of the initial dense torus at 100-R 9 ), and £ or b = 2-7r/f2o is the 
orbital time at Rq. The scale over which mass is deposited is 
10R g and the mass source function falls steeply beyond that. 
The mass is added with the local velocity and internal energy 
density. 

The initial condition of our RIAF to thin disk transition 
run (Fig. 5(a)) is the end state (at 24 orbits) of our low den- 
sity run with po = 10 9 m p g cm" 3 and a — 0.01 (similar to 
Fig. 3(a)). Figure 5(b) shows how the addition of mass looks 
like at early times — the mass being added at the local tem- 
perature (T ~ 10 8-9 K) in a circular region of radius 10R g 
about -Ro- Since mass addition is continuous, the density of 
the flow starts increasing, and after a few orbits the cooling 
time becomes shorter than the viscous time. This leads to 
the formation of a cold thin disk having T ~ 10 4 K. With 
time, as the density increases further because of mass addi- 
tion, there is more efficient cooling of the hot coronal gas into 
a geometrically thin disk which extends both inwards and 
outwards. The final steady state shows a geometrically thin 
disk extending very close to the black hole, surrounded by a 
hot low density corona with T ~ 10 9 K. The expansion of 
the thin disk towards the black hole takes place at the vis- 
cous timescale (tvisc) at the inner edge of the disk in Figure 
5(c). With more mass added, the disk expands viscously both 
inwards and outwards (Figs. 5(d,e)). 



3.3.1 RIAF to Thin Disk Transition 

A geometrically thin cold disc can only arise from a hot RIAF 
if the ratio t CO oi/ivisc < 1. Such a cooling-dominated dense 
state can be brought about by an increase in the mass ac- 
cretion rate. This can happen due to the variability of wind 
from the massive companion (e.g., Eversberg, Lepine, & Mof- 
fat || 1998 1, due to instabilities close to the Roche lobe (e.g., 
Sawada, Matsuda, fc Hachisu |1986 1 , or due to viscous insta- 
bility of accretion flows (say, due to thermal/viscous instabil- 
ity cause d by fast opacity variations with t emperature close 
to 10 4 K; Meyer fc Meyer-Hofmeister |l981 1. A recent exam- 
ple of the variability of mass accretion rate on large scales 
in AGNs is the discovery of a dense gas cloud falling in the 
accretion zone of Sgr A*, the Galactic center BH; this cloud 
is expected to cause a sudden increase in the mass accretion 
rate followed by a significant brightening of its X-ray and 
broadband emission (|Gillessen et al. 20121. Of course, the 



amount of mass dumped in this case is too small to cause the 
formation of a thin black-body disk. 

The mass and the density of RIAFs in our simulations 



3.3.2 Thin Disk to RIAF Transition 

The transition from a geometrically thin cold disk to a hot 
RIAF requires the reduction of mass in the accretion flow such 
that tcooi/ivisc J> 1 throughout. This happens naturally if the 
mass accretion rate in the thin disk at inner radii is larger than 
the mass addition rate at larger radii. Thus the thin disk to a 
RIAF transition happens on the viscous timescale at the outer 
thin disk (mass is dominated by outer radii), which can be 
quite long. Transient energetic events such as large magnetic 
flares and explosive nuclear burning are not expected to get 
rid of sufficient mass from the thin disk because most of the 
energy is expected to be coupled to the hot corona and not 
the cold massive disk. 

Figure 6 shows the time evolution of a high density run 
with a steady thin disk. The initial condition of this run is 
obtained by evolving our standard po = 10 11 7n p g cm -3 and 
a = 0.1 run for 24 orbits. Figure 6 shows the temperature 
contour plots as transition occurs from a geometrically thin 
to a geometrically thick accretion flow. This transition hap- 
pens simply because of exhaustion of mass due to accretion. 
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Figure 6. Transition from a geometrically thin to a geometrically thick flow due to viscous exhaustion of the thin disk. The grid resolution 
is 256 2 . 



The thin disk slowly recedes from the inner region until disap- 
pearing completely, giving way to a hot geometrically thick, 
low density corona (RIAF). The viscosity parameter a is cho- 
sen to be 0.1 instead of a — 0.01 in order to accelerate the 
evolution, which still takes about 40.5 orbits to transform into 
a RIAF. 



3.4 Density Upper Limit in Corona 

In the previous sections we have established the condition 
necessary for the formation of a cold, geometrically thin disk, 
namely, t C ooi/ivisc < 1. We have also observed that if the ini- 
tial density is high enough to trigger bremsstrahlung cooling, 
the final steady state consists of both a geometrically thin cold 
disk and an inner hot low density corona. The next question 
which arises is about the distribution of matter in the hot and 
cold phases. More precisely, what is the gas density distribu- 
tion in the hot corona around a cold thin disk? 

In order to address this question we have carried out 
simulations with different forms of accretion viscosity: v oc 
r 1 ' 2 , v oc r°, and v oc p. The magnitude of the viscosity is 



chosen such that t coo i/t v i sc <C 1 in the RIAF and hence a cool 
disk condenses. A subtle point to note is that this ratio should 
be calculated in the RIAF state before a thin disk condenses. 
The ratio t coo i/t v i sc can be ^C 1 in the corona at the radius 
where a thin disk exists because here the corona is in thermal 
equilibrium (e.g., Figure 7(b)). Any imbalance between the 
cooling and heating of the corona will be compensated by 
condensation on (or evaporation from) the thin disk. Figure 
7 compares the ratio of cooling to viscous time (t CO oi/ivisc) and 
the ratio of cooling to free-fall time (t coo i/iff) at late times, 
for these runs with different forms of viscosity. We observe 
that the ratio t C ooi/ivisc in the inner hot corona surrounding 
the cold thin disk is different for the different forms of v\ this 
is expected since t v i sc depends on v. However, the minimum 
value of tcooi/iff in the corona close to the thin disk, for all 
the three cases, is quite similar (~ 10 — 100). Thus we can 
conclude that the ratio of the cooling time to free-fall time, 
tcooi/ts ~ 10 — 100, corresponds to the upper limit on the 
density of the hot corona, irrespective of viscosity. 

The corona is cooling due to bremsstrahlung and is 
heated due to viscous dissipation (Eq. (3)). Thus it is in rough 
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thermal balance in steady state. Even if the corona is in ther- 
mal balance, mass cannot be added to it indefinitely, as a 
dense corona is expected to undergo multiphase cooling due 
to local thermal instability; excess mass should condense and 
fall onto the cold disk. This is analogous to what happens in 
dense cores of galaxy clusters, which, as shown by |Sharma et| 
al. | ( |2012[ ), self-adjust such that the ratio of the thermal in- 
stability timescale (£ti) to the free- fall timescale [tg) satisfies 
iTi/iff J> 10 everywhere, in order to achieve global thermal 
balance. The minimum t coo \/ts ratio in our accretion disk 
corona is 10 — 100 and respects the above limit, irrespective 
of viscosity. A detailed investigation of subtleties due to dif- 
ferences in the cooling time and the thermal instability time 
spherical versus cartesian gravity (see McCourt et al. ||2012 



Sharma et al. 20121, etc. are beyond the scope of this pa- 



per. The heated corona can lead to outflows but we have not 
evolved our simulations long enough to verify this. 

The distribution of gas in the disk vs. the corona is fur- 
ther illustrated by Figure 8(a), which shows the mass of the 
cold (T ~ 10 4 " 5 K) and hot (10 5 sC T < 10 105 K) gas in the 
accretion flow, as po is varied. Note that the cold gas repre- 
sents the geometrically thin disk, while the hot gas makes up 
the low density, geometrically thick corona. We observe that 
for po < 10 10 m p g cm -3 there is no cold gas. However, as po 



increases and crosses the t CO oi/ivisc < 1 threshold the mass of 
the cold disk increases steeply. On the contrary, the mass of 
hot gas increases much slowly compared to the mass of the 
cold gas. As we have already discussed, the coronal density 
has an upper limit (corresponding to t coo \/tg > 10 — 100), 
so the increase in the hot gas mass with increasing po occurs 
because of a bigger size of the corona, and because of a larger 
fraction of gas at higher densities (and lower temperatures). 

3.5 Transition Radius 

An important indicator of the structure of the accretion flow 
is the transition radius (Rtrans), the radius at which the hot 
accretion flow at smaller radii (R < -Rtrans) transitions to 
a thin disk at larger radii (R > -Rtrans). There is so far no 
accurate method to estimate the exact location of Rtrans, 



both theoretically and observationally (Narayan & McClin- 



tock 20081. For quantitative estimates we define Rtrans to 
be the radius where H/R ~ 0.02, recalling that the geomet- 
rically thick RIAF has H/R » 0.1 and H/R <; 10~ 4 for a 
thin disk. We have already seen that for our simulations with 
po < 10 10 m p g cm -3 there is no cold gas and the entire flow 
is a hot RIAF. (If the initial torus is placed much farther out, 
then a thin disk can form at larger radii even for a smaller 
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initial po-) Thus, the transition radius for po = W 9 m p g cm -3 
can be considered to be > 100R g . As the density of the flow 
increases (which is equivalent to an increase in the mass ac- 
cretion rate), we note that -Rtrans decreases monotonically, 
i.e., the cold thin disk moves closer to the black hole. Thus, 
a smaller .Rtrans signifies a cooling/thin-disk dominated flow, 
while a larger Rtrans indicates that the flow is dominated by 
a hot radiatively inefficient corona. 



4 A MODEL FOR BLACK HOLE TRANSIENTS 

In previous sections we established the condition for the for- 
mation of a thin disk and showed how cooling and viscous 
evolution brings about state transitions. In this section we 
show the results of a single simulation with a — 0.1 and ini- 
tial density po = 10 9 m p g cm -3 (Run A.ln9), which demon- 
strates the complete transition from a RIAF to a cold thin 
disk and back to a RIAF. Previous works, like ours, have 
identified the low-hard and high-soft states with a RIAF and 
a geometrically thin black body disk, respectively, but have 
not studied the time evolution of BH transients in detail (e.g., 
Esin, McClintock, fc Narayan |1997 |. 



We initialize a hot RIAF with T ~ 10 ' K everywhere 
(Figure 9(a)). This initial condition is obtained by evolving a 
po = 10 9 m p g cm~ 3 torus initial condition for several orbits. 
The cooling time is long in this case due to the low density, 
and hence there is no signature of a thin disk. Next, we add 
mass (after about 0.02 orbits at Ro) to the flow according to 



the density source term in Eq. 13 Figure 9(b) shows the onset 



of cold gas formation due to the increase in density and the 
decrease in cooling time brought about by mass addition. As 
the mass accumulates (much) faster than the accretion time, 
the density of the flow keeps increasing and eventually a thin 
cold disk (Figure 9(c)) is formed which extends very close to 
the black hole (<J 10R g ). The condensation of the initial thin 
disk occurs on the fast cooling timescale. The initial thin disk 
spreads inward on the viscous time at the inner edge of the 
thin disk in Figure 9(b). 

Once a steady thin disk is formed (after about 11 orbits), 
we stop the mass addition and the disk is simply allowed to 
evolve viscously. As accretion continues, sufficient mass is ex- 
hausted from the disk and the thin disk recedes far away 
from the central black hole (Figure 9(c)). As more and more 
mass is accreted into the black hole, the density of the flow 



© 2013 RAS, MNRAS 000,[T]{T5 



Black Hole Transients induced by Cooling 13 



50 



5£ 



mdot = 0.02 

hot 

mdot , = 0.09 

cold 



50 







mdot = 8.2e-4 

hot 

mdot ..= 5.5e-3 

cold 

-50' 



11 orbs 



(c) 



0.2 orbs 



* 



mdot = - 0.06 

hot 

mdot ,,,= 0.05 

cold 



(b) 



76 orbs 



(d) 



I 



t=0 



50 
R/R 



100 



(a) 



mdot = 6.6e-5 

hot 

mdot =0 

cold 



g 



50 
R/R 



100 



lo 9 10 T 



,12 



11 



10 



-9 



-7 



Figure 9. Complete state transition diagram, (a) Low-hard state (b) Intermediate hard state (c) High-soft state (d) Intermediate soft 
state. The intermediate soft state transitions to the hard state at > 100 orbs, m^ot and fn co \d are the accretion rates in units of M^dd for 
the hot and cold gas respectively. The grid resolution is 256 2 . 



drops, causing the cooling time to become longer than the 
viscous time, eventually resulting in the formation of a RIAF 
(in about 100 orbits from the start). The luminosity in this 
RIAF state continues decreasing on the viscous timescale at 
the circularization radius (~ Ro in the present case). Thus we 
return to the original state with which we started (similar to 
Figure 9(a)). Note that the thin disk state is the most long- 
lived of the transient states, since it takes a significant amount 
of time to remove sufficient mass from the disk through ac- 
cretion, as this process takes place at a viscous time at the 
outer edge of the steady thin disk, which is pretty long. Re- 
call that we are using v oc r ' 2 as the functional form for the 
viscosity, for both the thin disk and the RIAF. Whereas, in 
reality, the viscous diffusion coefficient for a thin disk is ex- 
pected to be (H/R) 2 ~ 10 -4 times smaller than in a RIAF. 
Thus, the black-body state is expected to be even longer lived 
than what is indicated in Figure 9. 

X-Ray astronomers popularly use the Hardness-Intensity 
diagram (HID) to describe the overall evolution of the spectral 
states of black hole XRBs, which are observationally known to 
be transient sources ( Remillard fc McClintock |2006 Fender 
fc Belloni |2 012). Typically in the HID, a black hole transient 
source transitions between various spectral states classified 



according to their X-ray luminosity and X-ray hardness trac- 
ing a 'Q'-shape, as shown, for example, by the 2002/2003 
outburst of GX 339-4 (|Belloni et al. ||2005| >. One can draw a 



one-to-one correspondence between the various spectral states 
and the ones shown in Figure 9 as follows — low-hard state 
(Figure 9(a)), hard intermediate state (Figure 9(b)), high-soft 
state (Figure 9(c)) and soft intermediate state (Figure 9(d)). 
This can be done, since Figure 9 also gives an estimate of 
the mass accretion rate (an indicator of the luminosity) cor- 
responding to each state. We must mention that the observed 
transients show short-timescale oscillations in hardness dur- 
ing the intermediate state which are probably not caused by 
variations in the mass accretion rate but are likely due to 
dynamical timescale magnetic processes in the corona (e.g., 
reconnection). 

We have already shown that the condition for the flow to 
be a RIAF is t cool /t viBC « 3nk B Tu/n 2 A(T)r 2 £ 1. This con- 
dition can be recast in terms of the accretion rate using the 
following scalings. From the equation of continuity one can 
write the number density n w M /2ni/m p H. The gas temper- 
ature for a RIAF is nearly virial and hence ksT ~ GMm p /r. 
Also, we write M in terms of the Eddington accretion rate 
Afedd = 4-KGMrrip/riaTC, where ot is the Thomson scat- 
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tering cross section and r\ = 0.1 is the radiative efficiency. 
Thus, the condition for the accretion flow to be a RIAF 
at 10R g , putting in all the numbers, is m < 0.1 a 2 , where 



m = M/MEdd (this was first pointed out by Rees et al. 
|1982[ ) . This indicates the existence of a critical accretion rate 
above which cooling will start dominating and a thin disk will 
form. 

The simulation in Figure 9 corresponds to a — 0.1. Thus, 
the condition for the flow to be a RIAF becomes m < 10 -3 , 
which is indeed the case in Figure 9. In all the panels rhcoid 
represents the accretion rate of the cold gas, while rhhot rep- 
resents the accretion rate of the hot gas. The 'accretion' rates 
are calculated by taking the difference of mass in the cold and 
hot phases at two consecutive snapshots. Figure 9(a) shows 
no cold gas, and rhhot = 6.6 x 10~ 5 <C 10 -3 ; this confirms 
that this flow is a low-luminosity RIAF with a low accre- 
tion rate (corresponding to the quiescent state). The state 
shown in Figure 9(b) corresponds to the hard intermediate 
state since mass in the hot gas becomes quite substantial, 
cold gas has just started to condense. The magnitude of the 
mass 'accretion' rate is w 0.001 (this is actually the mass ad- 
dition rate due to the source term in Eq 
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critical value for the formation of a thin disk (10 -3 
phase accretion rate rhhot is negative since the rate at which 
mass is added to the hot phase is smaller than the rate at 
which hot gas condenses into the cold phase. Figure 9(c) cor- 
responds to a steady thin disk state with rh co id = 0.09, which 
is greater than rhhot = 0.02; thus it is a cooling dominated 
flow. Also, rh « 0.1 3> 10 -3 indicates that this state has a 
higher luminosity and a prominent thin disk, namely it is the 
high-soft state. Figure 9(d) corresponds to the soft intermedi- 
ate state with rh co id = 5.5 x 10~ 3 , which is still greater than 
rhhot = 8.2 x 10~ 4 . While the thin disk has receded signifi- 
cantly in this state, there are still some signatures of the thin 
disk. 

Note that, although in our simulations we mainly study 
the geometrical structure of the accretion flow, we can still 
get an idea about the luminosity and the degree of hard- 
ness/softness of emitted X-rays. The soft photons are radiated 
by a geometrically thin, optically thick multi-color blackbody 
disk. The hot RIAF, on the other hand, is associated with 
hard X-rays, and the corresponding spectrum is dominated 
by a power-law; the power law arises because of the extension 
of the corona (and its temperature) over several decades in 
radius. Note that the conversion of the mass accretion rate 
(m) into a luminosity is not straightforward. The thin disk 
luminosity is ~ GMM /2Ri n where R ln is the inner radius of 
the cold disk; it not only depends on M but also on the inner 
truncation radius of the thin disk. Similarly, the emissivity of 
the RIAF is n 2 A(T), so the RIAF luminosity is proportional 
to m 2 and subdominant compared to the energy released via 
accretion ~ GMM/27?lso (where -Rlso is the last stable or- 
bit); most of the RIAF accretion energy goes into outflows 
(e.g., Fig. 1 in |Churazov et al. |[2005| . 

From the above considerations we note that Figure 9 de- 
scribes a complete state transition, (a) — >• (b) — > (c) — > (d) — > 
(a), which can be explained by variability in the mass accre- 



tion rate and viscous evolution of the accretion flow. Thus, 
very interestingly, our model for state transitions qualitatively 
reproduces the features of the HID. The observed timescales 
in lightcurves and hardness evolution can be used to constrain 
the circularization radius of the transient thin disk, its outer 
radius, the transition radius, and the total mass added in the 
accretion event. 



5 CONCLUSIONS 

The key conclusions of the paper are: 

• The condition for the transition of a RIAF to a cold, 
geometrically thin disk is that the cooling time be shorter 
than the viscous time, i.e., t C ooi/ivisc < 1 (see Fig. 3). When 
the density of the accretion flow exceeds the critical density 
corresponding to i CO oi/ivisc jj 1, cooling processes such as 
bremsstrahlung are triggered, resulting in a globally stable 
two-phase steady state consisting of an inner hot RIAF and an 
outer geometrically thin disk. The accretion flow with cooling 
is sensitive to the Shakura-Sunyaev viscosity parameter a. 
A higher a, which implies a shorter i v isc, requires a higher 
density to trigger cooling and hence to form a thin disk (see 
Fig. 4). 

• While whether a cool disk condenses is determined by 
icooi/ivisc, the ratio of the cooling time and the free-fall time 
determines the structure of the disk corona. The minimum 
ratio of the cooling time to the free-fall time in the corona 
close to the thin disk is found to be t coo \/tg ~ 10 — 100 (see 
Fig. 7). This ratio is independent of the kinematic viscosity u 
and indicates an upper limit for the density of the hot coronal 
phase, which is in thermal balance. If mass is added to the 
corona such that this density threshold is exceeded, the excess 
mass simply condenses and falls onto the thin disk. In the 
thin disk state the accretion flow has a subdominant but a 
substantial corona, which should produce some hard radiation 
even in the high-soft state (see Fig. 8(a)). 

• With the increase of the mass accretion rate, the transi- 
tion radius -Rtrans between the inner hot RIAF and the outer 
thin disk decreases monotonically, i.e., the thin disk moves 
closer to the black hole with an increasing mass accretion 
rate (see Fig. 8(b)). 

• We have qualitatively simulated transitions between vari- 
ous X-ray spectral states of black hole transients and their rel- 
ative positions in the Hardness-Intensity diagram/ Q-plot (see 
Fig. 9). Our model gives an estimate of the overall timescales 
in XRB transients. From our model we conclude that the tran- 
sition from a hot RIAF to a geometrically thin disk occurs 
due to an increase in the mass accretion rate. This transi- 
tion progresses on the viscous timescale at the inner edge R[ n 
of the initially formed thin disk (which can be quite fast if 
the circularization radius is small). The timescale of transi- 
tion from a low-hard to a high-soft state can thus be esti- 
mated to be ivisc ~ R 2 n /ac s H. For a 10M© BH with the 
inner disk at Ri n ~ 10 3 i? 9 , a — 0.1, the sound speed of 10 5 
K cold disk c s ~ 3 x 10 6 cm s" 1 and H/Ri n ~ 0.01, we es- 
timate t v isc,R in ~ 10 days. Similarly, our model shows that 
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the transition from a cold, thin disk to a geometrically thick, 
hot RIAF occurs at the viscous timescale at the outer edge 
-Rout of the steady state thin disk (which is much longer). 
For the same XRB considered above, if -Rout is ~ W 4 R g , we 
obtain t v isc,R out ~ 1 year (here we have assumed that tem- 
perature T oc r -3 ' 4 ). The increase in luminosity in the hard 
state ((a)— >(b) in Fig. 9) should happen on the mass addition 
timescale, which can be pretty short (of the order of dynami- 
cal time). These timescales match pretty well with generic ob- 
servations of BH XRBs. Since Roche-lobe filling XRBs have 
larger circularization radii (and hence longer viscous times) 
compared to the XRBs accreting massive stellar winds, the 
global timescales in the Q-plot for stellar wind accreting BH 
XRBs should be shorter. 
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